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Abstract 

Sudoku puzzles are now popular among people in many countries across the world with simple 
constraints that no repeated digits in each row, each column, or each block. In this paper, we demonstrate 
that the Sudoku configuration provides us a new alternative way of matrix element representation by using 
block-grid pair besides the conventional row-column pair. Moreover, we discover six more matrix element 
representations by using row-digit pair, digit-row pair, column-digit pair, digit-column pair, block-digit 
pair, and digit-block pair associated with a Sudoku matrix. These parametric Sudoku associated matrix 
element representations not only allow us to denote matrix elements in secret ways, but also provide us 
new parametric two-dimensional bijective mappings. We study these two-dimensional bijections in the 
problem of image scrambling and propose a simple but effective Sudoku Associated Image Scrambler only 
using Sudoku associated two dimensional bijections for image scrambling without bandwidth expansion. 
Our simulation results over a wide collection of image types and contents demonstrate the effectiveness 
and robustness of the proposed method. Scrambler performance analysis with comparisons to peer 
algorithms under various investigation methods, including human visual inspections, gray degree of 
scrambling, autocorrelation coefficient of adjacent pixels, and key space and key sensitivities, suggest 
that the proposed method outperforms or at least reaches state-of-the-art. Similar scrambling ideas are 
also applicable to other digital data forms such as audio and video. 
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I. Introduction 

In the new century, business news about conventional industries and companies keep reminding us how 
much impact the digital revolution bring to us and to the world. United States Postal Service is struggling 
hard to survive these days and have already closed hundreds of facilities and sites, because paper letters 
are not used as often as they were in the last century after email and other electronic communications 
dominates the way that people communicate with each other. Borders, an international book and music 
retailer once pioneered in the kind of large supermarket style bookstores, failed to survive under the 
impacts from online and e-books and has already bankrupted in the early of the year. Similar stories also 
happened to Eastman Kodak company, commonly known as Kodak, an American multinational imaging 
and photographic equipment, materials and service company that once dominated the photographic film 
industry for many years and invent and built the first digital camera as early as in 1975 fll, but now was 
bankrupted due to its slowness in transitioning to digital photography. As emails, digital images, digital 
books and other digital data carriers play an important role in people's life, the demand on the secure 
storage and transmission on digital data become a problem must be solved. 

Digital images, one typical type of two dimensional data, are considered to contain a huge amount 
of information, for example, a family photo might tell not only who are in this family and how they 
look like, but rough ages for each member and maybe their healthy conditions; a diagnosis CT image 
might tell a doctor whether the patient is healthy or not and if he/she is sick how bad he/she is; and 
a satellite photo might give information whether the interested region is under constructions and what 
these constructions could be. Because the information contained in a digital image and the information 
might be inferred beyond a digital image, it is very important to protect these information from any 
unauthorized use. One way of protecting digital images is called image scrambling, which disorders 
pixel relationship in the original image so that the scrambled image with rearranged pixels become 
unintelligent and unrecognizable. 

Mathematically speaking, all image scrambling algorithms rely on a two dimensional bijection, which 
maps pixels in the original image to rearranged pixels in the scrambled image. According to the generation 
of a two dimensional bijection used in a scrambler, we can classify a scrambler into groups: chaos based 
|[2), |[3|, spatial transform based ||4|, |[5|, matrix decomposition based |6| and cellular automata based Q, 
|[8|. According to the way that an image scrambler uses a two dimensional bijection, we can roughly 
classify it into two groups: conventional scrambling [j6J and total scrambling |[2|, (4), ||7|, |[8|. Conventional 
scrambling consider each pixel an unit and only shuffles pixel positions, while total scrambling consider 
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each bit of a pixel an unit and shuffles pixel bit positions and pixel positions. Generally speaking, the total 
scrambling scheme is considered more secure than the conventional scheme because image histograms 
have been changed after total scrambling, but stays unchanged after conventional scrambling. Moreover, 
image scrambling are an essential component for advanced data protection techniques, for example, image 
encryption |[9|-|[T8| and watermarking |[T9l-|[22|. 



In this paper, we extend our work in pT] , p3] | and introduce new ways of constructing two dimensional 
bijections for digital image scrambling using Sudoku matrices. We demonstrate that a N x N Sudoku 
matrix can be used to define six new parametric matrix element representations for a x matrix. 
Consequently, we are able to construct Sudoku associated two dimensional bijection by mapping one 
matrix element representation to the other. Moreover, we demonstrate that Sudoku associated two dimen- 
sional bijections actually performs scrambling in a guarantee way. For example, the Sudoku associated two 
dimensional bijection mapping from the row-digit pair to the row-column pair is to scramble image pixels 
in such a way that none two pixels originally lies in the same column will be in different column after 
scrambling, and the Sudoku associated two dimensional bijection mapping from the digit-column pair to 
the row-column pair is to scramble image pixels in such a way that none two pixels originally lies in the 
same row will be in different row after scrambling. In addition, we propose a total scrambling scheme 
Sudoku Associated Image Scrambler only using on these Sudoku associated two dimensional bijections 
for scrambling. Simulation results on various image types and datasets demonstrate the effectiveness 
and robustness of the proposed method. Visual and Analytical comparisons to peer algorithms suggest 
that the Sudoku Associated Image Scrambler outperforms or reach state-of-the-art of image scrambling. 
Statistical testing results also support that Sudoku Associated Image Scrambler successfully decorrelates 
pixels in the original image to random-like. 

The commonly used symbols and notations in this article are given in Table [l| In the rest of the 
article. Section II gives a brief background about Sudoku matrices; Section III discusses the Sudoku 
associated matrix element representations and two dimensional bijections; Section IV proposes the Sudoku 
Associated Image Scrambler using Sudoku associated two dimensional bijections; Section V shows 
extensive simulation results and compares the performance of Sudoku Associated Image Scrambler with 
peer algorithms in detail; Section VI concludes the article. 

II. Background 

The name Sudoku is the abbreviation of the Japanese 'Sunji wa dokushin ni kagiru\ which means 
'single number'. Conventionally, Sudoku refers to a number-based puzzle, consisting of 9 x 9 grids^ 
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TABLE I 

Symbol and Notation Usage 



Symbol 


Description 


N 


a squared integer 


H-) 


an indication function on range {0, 1} 


r 


row index of an element in a matrix 


c 


column index of an element in a matrix 


b 


block index of an element in a Sudoku matrix 


9 


grid index of an element in a Sudoku matrix 


d 


digit index of an element within a Sudoku matrix 


Xr,c 


a pixel/element located at the intersection of the rth row and cth column of image/matrix X 


s 


a Sudoku matrix 


I 


the set of natural numbers from 1 to 


/ 


a two dimensional bijection on I x I 


r' 


the inverse bijection of / 




a two dimensional bijection mapping from a representation R to R' 


Ir^r' 


a two dimensional bijection mapping from a representation R' to R 


fix{.) 


rounding function towards to zero 


rem(.) 


remainder function 


F 


a fixed matrix element representation pair 


P 


a matrix element parametric representation pair 


Rf 


a fixed matrix element representation using pair F 


R% 


a matrix element representation using pair P with reference Sudoku S 


o 


the function composition symbol 


W 


the width of an image 


H 


the height of an image 


T 


the number of pixels of an image 


K 


a scrambling key 


r(.) 


the Gamma function 


M 


the absolute value function 


9{-) 


the Student's t-distribution function 



divided into nine 3x3 blocks [ [24| . The objective is to complete the grids using digits ranging from 1 
to 9, in a manner that there are no repeated digits in any single row, column and block of the overall 
puzzle. 

Fig. [T] shows a Sudoku puzzle in a newspaper and its solution. The Sudoku hints given [T}(a) are 
identified by black colored numerals; the nine Sudoku block indices are identified by the large blue 
colored numerals ranging from 1 to 9; and the blank Sudoku elements in [T]-(a) are identified by red 
colored numerals in Fig. [T}(b), respectively. This is a conventional Sudoku puzzle, with a 9 x 9 size, to 
be filled with digits ranging from 1 to 9, and divided in square blocks of size 3x3. 

Throughout of this paper, we are interested in a Sudoku solution instead of a Sudoku puzzle, because 
multiple Sudoku puzzles can be made with respect to a single Sudoku solution, and it is the Sudoku 
solution that satisfies all digit constraints along rows, columns and blocks. Specifically, we consider a 
Sudoku solution from the point of view of matrix and thus call it a Sudoku matrix. It is noticeable that 
one can easily extend the conception of a 9 x 9 Sudoku matrix for other squared sizes, e.g. 4x4, 16 x 36, 

4n some literature, this 3x3 block is referred to as a box, or a square. 
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Fig. 1. A newspaper style Sudoku puzzle and its solution, (a) a Sudoku puzzle, and (b) its solution 
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Fig. 2. Large size Sudoku matrices, (a) a 25 x 25 Sudoku, and (b) a 36 x 36 Sudoku 



25 X 25 etc. Fig. [2] shows sample results for large size Sudoku matrices. Literally, we are therefore able 
to define a N x N Sudoku matrix with = as shown in Def. 1. 

Definition 1 (Sudoku Matrix): A N x N matrix is a Sudoku matrix if and only if its elements satisfy 
all three constraints that matrix elements in any row, and in any column, and in any nxn block contains 
exact N digits from 1 to A^. 

A N X N Sudoku matrix conforming Def. 1 has the properties including, but not limited to the three 
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listed below 

• The N elements within each row of this Sudoku matrix is a permutation of the nature number 
sequence {1, 2, • • • , N} 

• The N elements within each column of this Sudoku matrix is a permutation of the nature number 
sequence {1, 2, • • • , N} 

• The N elements within each block of this Sudoku matrix is a permutation of the nature number 
sequence {1, 2, • • • , N} 

• A subset of X Sudoku matrices can be perimetrically generated [ [TT| 

It is obvious that these three properties are directly from the definition of a Sudoku matrix. 

III. Methodology 

In this section, we mainly explore the possible matrix representations associated with a Sudoku matrix. 
And we show that element-wise mapping between any of two representations is actually a bijection, which 
ensures to shuffle and restore image pixels in an easy way. 

A. A Sudoku Matrix Definition 

Mathematically, we can define a N x N squared size Sudoku matrix via an indicator function / as 
follows: 

f 1, if Xrc = d 

Ix{r,c,d)=l ' (1) 

I 0, Otherwise 

where r, c and d denote the row index, column index, and digit index of a matrix element in matrix X; and 
Xr^c denotes the element at the intersection of the rth row and cth column in X. Let I = {1, 2, • • • , N} 
be the index set. Then 

Definition 2 (Sudoku Matrix): A N x N matrix 5 is a Sudoku matrix if and only if 

N 

• For arbitrary r, d G I, we have Yl Is{r^ c, rf) = 1 

r=l 
N 

• For arbitrary c, (i G I, we have Yl ^si^^ c, d) = 1 

c=l 
N 

• For arbitrary 6, d G I, we have ^ Isi^r^ bc^d) = 1 

where r, c, b and d denotes the row index, column index, block index and digit index, respectively; br 
and be denote the row index and the column index of the ^th grid in the 6th block. 
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In particular, the row index r and the column index c of the ^th grid in the 6th block are defined as 
follows, 

{b^ — rem{h — 1, n) • n + rem{g — 1, n) + 1 
(2) 
be = fix{b — 1, n) • n + fix{g — 1, n) + 1 

with the rounding function towards to zero fix{.), namely fix{p, q) = [p/q\ and the remainder function 
rem{.) , namely rem{p^ q) = p — fix{p^ q) • q. In this way, we define the two-dimensional mapping from 
the Sudoku representation R[b,g) to the conventional representation R(^r,c)- It is noticeable that this is a 
bijective mapping, and we can write the inverse mapping as 

{g — rem(r — 1, n) • n + rem{c — 1, n) + 1 
(3) 
b = fix{r — 1, n) • n + fix{c — 1, n) + 1 

Although mathematically the bijective mapping between representation R(^r,c) representation R{b,g) 
can be explicitly written as Eqs. ([2]) and ([3]), such bijective mapping is actually nothing but to link two 
representations with respect to each matrix element. 

Fig. [3] shows an example of matrix element representations R(^r,c) the R{b,g) for a 4 x 4 matrix. As 
can be seen, the actual mapping from representation c) to representation R(b,g) defined in Eq. Q can 
be simply described to map a (r, c) pair in Fig. [3]-(b) to its corresponding (6, g) pair in Fig. [3]-(c) with the 
same color. Similarly, the actual mapping from representation R{b,g) defined in Eq. ([3]) to representation 
R(^r,c) is to map a (6, g) pair in Fig. [3]-(c) to its corresponding (r, c) pair in Fig. |3]-(b) with the same 
color. Consequently, the bijection between these two representations are shown in Fig. |3]-(d). 

So far we have two matrix element representations, namely row-column pair representation R(^r,c) 
block-element pair representation R(b,g)' To differentiate from those representation associated with one 
or more Sudoku matrices, we call these two representations fixed representations because they are fixed 
rather than parametric and use symbol F to denote a fixed representation pair, i.e. F G {(r, c), {b^g)}. 

B. Sudoku Associated Matrix Element Representations 

Considering the task of digital image scrambling, the two dimensional bijection fn^^^^^^R^^^^g^ and 
fji^^ are inappropriate in the sense that such bijection is fixed. With the help of a Sudoku matrix, 

however, we are able to easily define two dimensional parametric bijections between two matrix element 
representations. 

Due to the digit constraints in each row, each column and each block within a N x N Sudoku matrix 
S, we are able to create parametric matrix element representations for a x matrix using the reference 
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Fig. 3. Matrix element representation example; (a) a 4 x 4 matrix with each element identified in a unique color; (b) the 
conventional row-column representation R{r,c)\ (c) the Sudoku representation R{h,g)\ and (d) the bijective mapping between 
matrix element representation i?(r,c) and R(b,g)- 



Sudoku matrix S. Specifically, we are able to have six parametric matrix element representations using 
row-digit pair (r, rf), digit-row pair (d, r), column-digit pair (c, d), digit-column pair (d, c), block-digit 
pair (6, d) and digit-block pair (d, b). Intuitively, one can convince him/herself that matrix elements can 
be represented with the (r, d) pair, because each matrix can be decomposed with respect to rows, and 
each element in a row is with a unique digit associated with a reference Sudoku in its row, implying 
that one can unique locate a matrix element whenever he/she knows the row index r, the digit index d 
and the reference Sudoku matrix S. Similarly, it is not difficult to see that other listed Sudoku associated 
pairs are also valid matrix element representations. We denote these Sudoku associated matrix element 
representations as Rp where S is the associated Sudoku matrix and P denotes a Sudoku associated 
representation pair, i.e. P G {(r, d), (d, r), (c, d), (d, c), (6, d), (d, b)}. 

Fig. |4] shows the six Sudoku associated matrix element representations for a 4 x 4 matrix, when the 
reference Sudoku S shown in Fig. |4]-(b) is used. It is noticeable that each representation is able to uniquely 
locate an element in M. Details about why these Sudoku associated matrix element representations are 
valid are given in Appendix. 

Besides the existence of the six Sudoku associated matrix element representations, we also want to 
emphasize that for the same Sudoku associated matrix element representation, it could be very different 
when the reference Sudoku matrices are different. Here is a simple example illustrating this fact for 
interested elements in a 4 x 4 matrix. Fig. [5]-(a) shows our interested elements, and -(b), (c), (d) give 
three reference Sudoku matrices. Then for each interested element, we are then able to represent it with 
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Fig. 4. Sudoku associated matrix element representations; (a) a 4 x 4 matrix M with elements identified in distinctive colors;(b) 
reference Sudoku matrix S; (c) Rfr,d) representation; (d) Rfd,r) representation; (e) Rfc,d) representation; (f) Rfd,c) representation; 
(g) Rft,d) representation; and (h) Rfd,b) representation. 

various Sudoku associated matrix element representations. For example, the green element located at (2, 3) 
in the conventional row-column representation is denoted as (2, 2) by using the (r, d) representation pair 
associated with Sudoku Si, because the element located at row index r = 2 and digit index d = 2 
in Si is the green element. In the same manner, this green element can be also denoted as (4, 3) by 
using the representation R(d,b) associated with Sudoku ^3. And Table shows the six matrix element 
representations for interested elements associated reference Sudoku matrix S^i, 5^2 and ^3. As can be seen 
from this table, a matrix element can be represented with Sudoku associated representation pairs including 
(r, (i), (d, r), (c, d), (d, c), (6, d) and (d, 6); and a matrix element can be represented differently by one 
Sudoku associated matrix element representation, if its associated reference Sudoku matrix changes. The 
reason behind such differences is that different Sudoku matrices have different digits in a grid. Moreover, 
due to the fact that these Sudoku associated matrix element representations are sensitive to the reference 
Sudoku matrix, we are therefore able to construct parametric two dimensional bijective functions for 
image scrambling. 

It is clear that for each Sudoku matrix, there are six associated matrix element representations. Since 
any two matrix element representation denote the same set of matrix elements, a bijective mapping 
then can be constructed correspondingly like we showed in Fig. [3]-(d). Since we have two fixed matrix 
element representations, namely row-column pair and block-grid pair, we are able to construct two new 
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Fig. 5. A 4 X 4 matrix with colored elements and two Sudoku matrices; (a) interested matrix elements identified by colors; 
(b) Sudoku matrix 5*1 ; (c) Sudoku matrix S2 and (d) Sudoku matrix 5*3 • 

TABLE II 

Matrix Element Representations Associated with a Reference Sudoku Matrix 
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bijections denoted as /rs^ji^ by mapping from the Sudoku associated matrix element representation 
Rp to the two fixed representations Rf, and another two new bijections by mapping from the two fixed 
representations to the Sudoku associated matrix element representation denoted as /rs^r^, where S 
denotes the reference Sudoku matrix, R denotes a matrix element representation, P denotes a Sudoku 
associated representation pair with P G {(r, d), (d, r), (c, d), (d, c), (6, d), (d, 6)}, and F denotes a fixed 
matrix element representation with F G {(r, c), (b^g)}. In summary, the total number Sudoku associated 
bijections can be directly constructed between a Sudoku associated matrix element representation to a 
fixed element representation is 6 x 2 x 2 = 24. 

Fig. [6] shows a naive one step image scrambling using the twenty-four Sudoku associated bijections. As 
can be seen from these results, different Sudoku associated bijections scramble the original Trees image 
differently. For example, the bijections /rs and /rs shuffle image pixels within each 

row; the bijections /rs .0 and /rs .0 shuffle image pixels within each column; the bijections 
fps .0^^ ^ and fps .0^^ ^ shuffle image pixels within each 16 x 16 block. It is also noticeable that 
/rs ^^^^^^^^, /rs ^^^^^^^^, and /rs ^^^^^^^^ are the three most powerful bijections that scramble the 
original Trees image to almost random-like from the point view of human vision inspection. 
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C. More Sudoku Associated Bijections 

It is well known that the composition of two bijective mapping is still a bijective mapping. In other 
words, if both /o and /i are two bijection on I x I, then fnew = /o o /i is also a bijection on I x I and 
so is fnew = /i o /o • Therefore, we can even find more two dimensional bijective mappings by simply 
composing two or more existent bijective mappings. 
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TABLE III 

The number of two dimensional BIJECTIVE mappings associated with SUDOKU MATRICES 



# of two dimensional bijection: Af^ 


# of Sudoku matrices: s 
12 m 


Composition times t 


t = 


24 


48 •• 


24m 


t = 1 


242 


482 .. 


(24m)2 


t = j 


24^'+i 


48^+1 . . 


(24m)^+i 



Table |nT] shows the relations of the number of reference Sudoku matrices and the number of two 
dimensional bijective mappings. When one Sudoku matrix is used for reference, then six new matrix 
element representations can be found. Therefore, when m Sudoku matrices are used, there are 6m 
Sudoku associated matrix element representations. Since a bijection can be constructed between a Sudoku 
associated matrix element representation and one fixed representation, the number of bijections available 
is 2 X 2 X 6m = 24m, where 6m is the number of Sudoku associated matrix element representations, 2 
is the number of fixe matrix element representations, and 2 implies the two possible ways to construct a 
bijection either mapping from the fixed representation to the Sudoku associated one or from the Sudoku 
associated one to the fixed one. 

Because the composition of two bijections is just a new bijection, we can define more bijections 
by involving function compositions. Without loss of generality, say we are interested in the number of 
bijections defined by t times of function compositions from i bijections /i, /2, • ' ' 1 fi^ 

fiooii---oit ~ fio ° fii ' ' ' ^ fit 

^ V ' 

t times 

where zq, ii, • • • , it G {1, 2, • • • , i}. Then it is not difficult to see that the total number of bijections by 
t times of compositions is since for each function fi. with j G {0, 1, • • • ,t} we have i candidate 
function and in total we have t -\- 1 functions to be determined. In Table [lll| we have 24m of possible 
Sudoku associated matrix element representations and j times of compositions i.e. i = 24m and t = j, 
and thus the number of bijections involving j times of compositions is i^^^ = (24m)^^^. 

It is worthwhile to note the number of bijections calculated in Table [lll| includes self-mapping, i.e. 
fnew = /o o /i with fi = /q~^. However, as long as the number of Sudoku matrices increases the 
possibility of randomly composing a pair of bijections (one is the inverse of the other) quickly approach 
to zero. An alternative remedy to avoid self-mapping is to choose the two Sudoku associated bijections 
with different reference Sudoku matrices, so that a new composed bijection from these two bijections is 
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impossible to be a self-mapping unless the two reference Sudoku matrices are identical. 

IV. Sudoku Associated Image Scrambling Scheme 

In previous sections, we showed that given N x N Sudoku matrix S for reference, there are six 
Sudoku associated matrix element representations, namely the row-digit pair representation R^^^y digit- 
row pair representation R^^^y column-digit pair representation R^^^y digit-column pair representation 
i?^ , block-digit pair representation i?^ and digit-block pair representation i?^ . Each Sudoku S 
associated matrix element representation Rp with P G {(r, d), (d, r), (c, d), (d, c), (6, d), (d, b)} and a 
fixed representation Rjr with F G {(r, c), (6,^)} can be used to construct a pair of bijections /rs^^r^ 
and fRS^^Rp, where Jrs^^r^ denotes the bijection mapping from the fixed representation Rp to the 
representation Rp, and /rs^r^ denotes the bijection mapping from the Sudoku associated representation 
Rp to the fixed representation Rp. Therefore, given a Sudoku matrix S and a fixed representation Rp, 
there are 24 Sudoku associated bijective mappings of /r^p and /r^jp- In this section, we use these two 
dimensional Sudoku associated bijections for digital image scrambling. 

A. Sudoku Associated Image Scrambler 

Although a Sudoku associated image scrambler can be designed in various means, we design a Sudoku 
associated image scrambler as shown in Fig. |7j where Key is an encryption key of 192-bit length. 
Parameter Generator generates key dependent and round dependent parameters including a Sudoku 
associated two dimensional bijection /rs^r^ or /rs^r^, Blockwise Scrambling shuffles pixels within 
3. N X N image block for every image block of the input image. Consequently, the image descrambling 
process is just the reverse of the scrambling process as shown in Fig. [7]-(b). In the following discussion, 
we consider an input image X is of size W x H with nB bit-depth. 

Parameter Generator translates a scrambling key K to required parameters for image scrambling. 
Specifically, it takes a scrambling key K to generate a set of parameters including twelve Sudoku 
associated matrix element representation pairs = {P^, ' ' ' 5 ^12}' twelve fixed representation pair 
= ' • • 5 ^12} the twelve mapping parameters = {m^^, m^r ' ' 1 ^12} the zth bit- 

plane in X as shown in Algorithm 1. In this way. Parameter Generator guarantees that for each bit-plane 
in X all six possible Sudoku associated matrix element representation pairs and two possible fixed 
representation pairs are used in the twelve round scrambling. Meanwhile, these parameters are bit-plane 
dependent, implying different bit-planes will be scrambled in a different way in the future processing. 
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Fig. 7. Flowchart of a Sudoku associated image scrambler; (a) scrambling process; (b) descrambling process 



Furthermore, Round Sudoku Generator takes the image size and the scrambler round u to generate a 
N X N Sudoku matrix as shown in Algorithm 2. 

Consequently, a two dimensional Sudoku associated bijection fl can be found for the ith bit-plane in 
uth scrambler round as follows 

" (4) 

[ J'Ki^f:^^ if = 1 

Once the X parametric Sudoku associated two dimensional bijection is obtained for the ith 
bit-plane in the uth scrambler round, we are then able to scramble image pixels with in x image 
within thoWxH input image for all of its bit-planes. In order to obtain pixel shuffling between blocks, 
we apply a image shift function imgShift{.) for each cipher round as shown in Eq. ([5]). 



Y — imgShift(X^ h) 



X{w ■.W,h:H) 


X{w ■.W,l:h-1) 


X{l:w-l,h:H) 


X{l:w-l,l:h-l) 



(5) 



Pseudo-code of Blockwise Scarmbling algorithm is given below in the MATLAB fashion. Consequently, 
we then shuffled every N x N image blocks within the input image X. 

In order to scramble an arbitrary W x H image without bandwidth expansion, we scramble blocks 
on the edge with a certain amount of overlapping with previous blocks as shown in Fig. [8] The block 
indices listed in Fig. [8] indicate the processing order of Blockwise Scrambling, implying to first shuffle 
image pixels in the upper-left block marked as index 1 in the figure, then block 2, block 3, so on and so 
forth until block 6. Since W and H may not be necessarily to be divisible by A^, we process blocks on 
the image edge with overlapping, e.g. block 2 and 3 are overlapped, and block 3 and block 6 are also 
overlapped. In summary, there are Oh x Ow pixels belongs to more than one block to be processed. In 



A PAPER DRAFT SUBMITTED TO IEEE TRANSACTIONS ON MULTIMEDIA 



15 



Algorithm 1. Parameter Generator (R', F', m') = BPG{K, i) 

Require: X is a key of 192 bits composed of twenty-four subkeys {K^^\ K^'^\ • • • , K^"^^^} 
Require: i is the ith bit-plane to be processed 

Ensure: (P, F, m) is a list of twelve pairs of Sudoku associated matrix element representation pairs, fixed matrix element representation pairs, 
and mapping directions. 

. . . , k(24)] = k % Divide key K to twenty-four subkeys K(^\ K(^\ • • • , ^(24) 
Q = [K*^ Krnod(i+i,24)+i^ . . . ^ ^mod(z+i2,24)+i] % Extract a twelve element subkey sequence starting at K^^^ 
idx = sort(Q) % Sort this sequence Q and generate the element index sequence in the sorted sequence 
for w = 1 : 1 : 12 do 

q = idx(j) % For each index in idx, pair it to a (R, F) representations 
if mod{q, 2) 1 then 
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end for 



Algorithm 2. Round Sudoku Generator Su = RSG{K, u, W, H) 

Require: K is a key of 256 bits 
Require: u is the cipher round number 
Ensure: Su is sl N x N Sudoku matrix 

Ku = R oundKeyG enerator{K, u) % Generate a round key ^ 

N = l^mm(W, if) J % Determine Sudoku size 

Su = SudokuGenerator{Ku, N) % Generate a N x N Sudoku matrix 
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Algorithm 3. Blockwise Scrambling Y — BlkScarmbling{X, fu, N, nB) 
Require: / is an image of size W x H 

Require: fu = {fu ^fu^''' ? fS^ } is a vector of nB bijections with each bijection for one bit-plane 
Require: N is a. squared integer 

Require: nB is the number of bit-planes contained in image X 
Ensure: Y is an scrambled image of X at the same size 

nRow = [W/N] 

nCol = \H/N] 

w = \W/12] 

h = \H/12] 

X = imgShift{X, w, h) % shift image X w pixels along rows and h pixels along columns 
for z = 1 : 1 : nRow do 
\f i = nRow then 

rO= {i — 1)N -\- 1 % 1st row index of the current block 
rl = iN % last row index of the current block 
else 

rO = W — N -\- 1 % 1st row index of the current block 
rl = W % last row index of the current block 
end if 

for j = 1 : 1 : nCol do 
if j = nCol then 

cO = {j — 1)N -\- 1 % 1st column index of the current block 

cl = jN % last column index of the current block 
else 

cO = H — N -\- 1 % 1st column index of the current block 
cl = H % last column index of the current block 
end if 

tBlk = I{rO : rl, cO : cl) % extract the current block 

tmp = % create a temporary intermediate variable 

for / = 1 : 1 : do 

tBP = bitget{tBlk, I) % extract the uth bit-plane 

oBP(rO : rl,cO : cl) = fl{tBP) % output scrambled bit-plane 

tmp = tmp + oBP •2}~^ % store results in temporary variable 

end for 

Y{rO : rl, cO : cl) = tmp % output scrambled block 
end for 
end for 



this way, we are able to fit a two dimensional Sudoku associated bijection defined on I x I for an input 
image with an arbitrary size W x H without bandwidth expansion. 

In regarding of descrambling, we can simply substitute a certain Sudoku associated bijection with 
its inverse. In other words, if the bijection used for scrambling is fRS^R^,, then the bijection used for 
descrambling is fRS^R^, and vice versa. 

B. Discussions 

Sudoku Associated Image Scrambler is a multiround scrambler that shuffles image pixels with respect 
to one Sudoku associated bijection for each round. This multiround scrambling process is equivalent to 
create a new bijection by taking the function composition of multiple bijections, which is discussed in 
Section III-C. It is not difficult to validate such equivalence between the multiround Sudoku associated 
image scrambler and the bijection by composing multiple bijections. Without loss of generality, denote the 
twelve two dimensional Sudoku associated bijections for the zth bit-plane are /I, /2, • * * ? fh Sudoku 
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Fig. 8. Overlapped blocks in Blockwise Scrambling 
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Associated Image Scrambler. Then running the scrambler for multiple times is nothing but to construct 
a new bijection fmw as follows 



fnen,{x, y) = fl,{fi,{- ■ ■ mi{x, y)))) = o /{^ • • • o o fiix, y) 



(6) 



which is actually the composition of the twelve bijections. 

Meanwhile, we want to emphasize that the Sudoku associated bijections used in Sudoku Associated 
Image Scrambler from one round to another is not independent, but dependent in the sense that in the 
twelve round the fixed representation R{r,c) -^(b,g) appears exact six times respectively and the six 
Sudoku associated representations appears exact twice for each. We set up this constraint to guarantee the 
strong scrambling effect. Although different Sudoku associate bijections have different scrambling impact 
on an image, many of them are of limited scrambling effect in the sense that scrambled images are not 
random-like as shown in Fig. |6j Therefore, if we allow random Sudoku associated bijections for each 
round, it is possible to pick twelve identical Sudoku associated matrix element representation pairs, i.e. 
Pi = P2 = • • • = P12 and also twelve identical fixed representation pairs = F2 • • • = ^12- Although 
the reference Sudoku changes for each round, the scrambled image might not change that much. Fig. [9] 
shows an example of using bijection /^s ^^^r^^ for image scrambling. As this mapping implies, it only 
scrambles pixels within a Sudoku block. And thus scrambled image gets blurry soon after applying this 
bijection once. Consequently, the scrambled image after applying this bijection twelve times shown in 
Fig. |9]-(c) is more or less the same as the previous, indicating that simply composing a number of Sudoku 
bijections does not necessarily improves scrambling quality. Fig. |9] -(d) shows the scrambled image with 
additional image shifting operation defined in Eq. (|5]). This result shows that this operation does help 
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improve scrambling quality by generating a more evenly distributed scrambled image. Fig. [9] -(e) shows 
the scrambled image using the twelve round Sudoku Associated Image Scrambler (in this particular case 
we consider image Trees to be bit-depth 1), which is completely unrecognizable and random-like. 




(a) (b) (c) (d) (e) 



Fig. 9. Image scrambling results discussion; (a) image Trees\ (b) scrambled image after one round /^s ; (c) scrambled 

image after twelve rounds fps ; (d) scrambled image after twelve rounds f^s with image shifting; and (e) 

scrambled image after twelve rounds scrambling in Sudoku Associated Image Scrambler. 

In regarding of scrambler security, we suggest to change a scrambler key K frequently. Because an 
image scrambler itself is not semantically secure and is vulnerable under certain attacks, for example, a 
simple but effective algorithm based on the chosen plaintext attack (CPA) model could be described as 
follows, 

1) Construct a all-zero image X of size W x H except the pixel Xij ^ 

2) Scramble this image X by using an image scrambler and a unknown key K and obtain scrambled 
image Y 

3) Find out nonzero value pixel in Y and record its position (vidx^Cidx) as the output of a bijection 
/, i.e. f{ij) = {ridx, cidx). where idx = (i - l)H + j 

4) Repeat the above three steps by changing the position of the none-zero pixel in X, for all i G 
{1,2,--- and jG{l,2,--- 

Consequently, the bijective mapping / is the equivalent bijection of the used scrambler for a x 
image under the key K. In other words, all x i7 images scrambled by the used image scrambler under 
key K can be perfectly cracked by inverse mapping scrambled image pixels using mapping f~^. And 
the complexity of the above attack isWxH. 

It is worthwhile to note that the above attack is a generic attack for all image scramblers. However, one 
remedy to this type of attacks is to change encryption key frequently for a scrambler. For example, say we 
use the Sudoku Associated Image Scrambler for high-definition television (HDTV) with the format 1080p, 
namely 1920 x 1080 pixels per frame. Then an adversary need to take 1080 x 1920 = 2073600 frames 
for fully crack an scrambling key K. We therefore set up a key change for every 518400 = 2073600/4 
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frames, which is equivalent to 6 hours of HDTV programs if the frame rate is 24 frames per second. As 
a result, an adversary may at most recover 25% of pixels in a HDTV frame, which is definitely of poor 
visual quality. 

V. Analysis and Comparison Results 

Digital image scrambling is to rearrange image pixels in a deterministic way but with a random-like 
appearance. In this section we focus our effects on performance analysis and comparisons for the Sudoku 
Associated Image Scrambler. 

A. Experiment Settings 

The following simulations are all performed under the Windows 7 operation system with the Intel Corel 
CPU @2.66GHz and QGB memory. In order to make easy comparisons, we run our Sudoku Associated 
Image Scrambler on test images (8-bit grayscale) whose scrambling results are widely reported and 
compare our results with other peer algorithms including the scrambling method based on two dimensional 
cellular automata proposed by Abu Dalhoum et al |8|, the scrambling method based on chaos map 
proposed by Ye |25|, the scrambling method of using cellular automata proposed by Ye et al ||7|, and 
the scrambling method based on ASCII code of matrix element proposed by Ye et al These test 
images includes Cameraman, Barbara, Lenna, and test images #157055, #69015 and #239096 in the 
Berkeley image segmentation dataset^ [ |26| . Since the results about these test images have been previous 
reported, we make comparison of the proposed scrambler with peer methods by analyzing these test 
images. 

Additionally, we also test the performance of the proposed scrambler on other image types than (8-bit 
grayscale) to validate that the proposed scrambler is able to deal with an image with an arbitrary bit 
depth, including binary images CCITT-3 and CCITT-7 from the CCITT fax image compression dataset 
16-bit grayscale knee MRI images, and color images 4.2.03 and 4.2.07 from the USC-SIPI miscellaneous 
image dataset ^. 

Details about these test images and results with analysis and comparisons are presented in next sections. 



^This dataset is free for non-commercial research and educational purposes, available under 



[Resea rch/Projects/CS/vision/bsds/ as the date of 04/10/2012. 

"^The CCITT database can be found under page: |http ://cdb.paradice- insight.us | as the date of 04/10/2012. 



http ://www.eecs .berkeley.edu/ 



The USC-SIPI image database can be found on 



http://sipi.usc.edu/database/ as the date of 04/10/2012. 
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B. Simulation Results 



Fig. 10 shows the simulation results of the proposed image scrambling method using the two dimen- 



sional Sudoku associated bijections with peer algorithms |[8|, |[25|. As can bee seen from these results, 



the proposed method outperforms the two recent algorithms, in the sense that its scrambled image is 
random-like and pixels more evenly scrambled, while results of Van De Ville et al. 's method |[6| contain 
distinguishable patterns of foreground from background; results of Abu Dalhoum et al. 's method ||8| 



contain line-like patterns and results of Ye's method p5| are of mesh-like patterns. Furthermore, it is 
noticeable that Abu Dalhoum et al. 's method might require extra space to deal with the edge effect in 
cellular automata (see results on the 2nd column), while the chaos map used Ye's method might be not 
that dynamic and generates weak results (see Ye's result on test image #239096). 



Besides the 8-bit grayscale images listed in Fig. 10 we also test the proposed method on various image 



types including binary images CCITT-3 and CCITT-7 from the CCITT fax image compression dataset ^, 
16-bit grayscale knee MRI images, and color images 4.2.03 and 4.2.07 from the USC-SIPI miscellaneous 
image dataset ^. Scrambling results are given in Fig. 



11 



In regard of the detail information about these test images including processing time and speed, we list 
these relevant information in Table |IV| Roughly speaking, the speed of the proposed scrambler is about 
150KB /s, or 0.15MB /s equivalently under MATLAB environment. It is worthwhile to note that this 
speed can be largely enhanced by using parallel computations, because the scrambling process of each 
bit-plane is independent of each other and image blocks that not overlapped are also independent of each 
other, both of which implying these works can be done in parallel to save time. Meanwhile, it is well- 
known that MATLAB is very slow for far-loop execution, and thus implementing the proposed scrambler 
in other languages might further enhance the processing speed. Therefore, the proposed scrambler with 
a proper implement could meet the demand for real time processings. 

C. Gray Difference and Gray degree of Scrambling 

Gray difference and gray degree of scrambling (GDD) |7| are two measures used for quality evaluation 
of scrambled images. The gray difference for aW x H image X is defined as follows, 

^^i^j ^4 X] ^^^^^ ~ ^i+p^j+qf ('7) 

MG{-1,+1} 
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Fig. 10. Image scrambling results with visual comparisons; 1st column: test images (from 1st row to the last are Cameraman, 
Barbara, Lenna, and #157055, #69015 and #239096); 2nd column: simulation results of Van De Ville et al. 's method |^; 
3rd column: simulation results of Abu Dalhoum et al 's method |^; 4th column: simulation results of Ye's method |25); and 
5th column: simulation results of ours. 
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Fig. 11. Image scrambling results on various image types using Sudoku associated image scrambler; (a) and (b) CCITT-3 
before and after scrambling; (c) and (d) CCITT-7 before and after scrambling; (e) and (f) knee MRI sample A before and after 
scrambling; (g) and (h) knee MRI sample B before and after scrambling; (i) and (j) 4.2.03 before and after scrambling, and (k) 
and (1) 4.2.07 before and after scrambling. 



Further, the mean gray difference of image X is computed by averaging all pixels except those on image 
edges as shown in Eq. ([9]). 

W-lH-l 

^^^^^ ^ iW-2){H - 2) 

Finally, the gray degree of scrambling is defines as the ratio of the difference of the mean gray differences 
before and after scrambling image X to their sum, namely 

Table |V| compares the gray degree of scrambling of the proposed method with peer algorithms ||4|, 
|[6|-|[8j, [[25t'^. As can be seen, the proposed method outperforms the listed peer algorithm by achieving 
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TABLE IV 

TEST Images Information with Execution Speed 



Image File 


Test Image Information 
Width Height Bit Depth 


Execution Information 
Time(s) Speed(KB/s) 


Cameraman 


256 


256 


8 


2.8280 


181.0467 


Barbara 


256 


256 


8 


2.9266 


174.9470 


Lenna 


256 


256 


8 


3.0336 


168.7764 


#157055 


321 


481 


8 


5.0176 


240.4053 


#69015 


481 


321 


8 


4.6216 


261.0044 


#239096 


321 


481 


8 


4.6537 


259.2040 


CCITT-3 


2339 


1728 


1 


31.4804 


125.3816 


CCITT-7 


2339 


1728 


1 


32.8107 


120.2980 


Knee MRI Sample A 


448 


448 


16 


21.1253 


148.4476 


Knee MRI Sample B 


448 


448 


16 


21.1838 


148.0377 


4.2.03 


512 


512 


24 


31.1277 


197.3805 


4.2.07 


512 


512 


24 


32.4009 


189.6244 



a higher gray degree of scrambling. 



TABLE V 

Gray Scrambling Degree of Scrambled Images 









Test Ima 


^es 






Method 


Cameraman 


Barbara 


Lenna 


#157055 


#69015 


#239096 


Van De Ville et al. 's 1 6 |, 2004 


0.1315 


0.1930 


0.1353 


0.2752 


0.5397 


0.4502 


Ye etal. 's |4|, 2007 


0.8832 


N\A 


0.9010 


0.8780 


0.8646 


0.8976 


Ye et al. 's |7 1, 2008 


0.8926 


0.8740 


0.9311 


0.8731 


0.8789 


0.9134 


Abu Dalhoum et al. 's 2012 


0.8971 


0.8749 


0.9320 


0.8821 


0.8827 


0.9388 


Ye's | 25 1, 2010 


0.9118 


0.9176 


0.9618 


0.9218 


0.9406 


0.9670 


Ours 


0.9165 


0.9200 


0.9663 


0.9359 


0.9439 


0.9672 



D. Adjacent Pixel Autocorrelation 

A digital image is commonly of high information redundancy in the sense that adjacent pixels are 
strongly correlated. In contrast, a well scrambled image should break such correlations between neighbor 
pixels, so that the scrambled image is unrecognizable. For example, one can simply plot a pixel correlation 
figure as shown in Fig.[T2j where each subfigure illustrating the adjacent pixel correlations before and after 
scrambling image Lenna by plotting the intensity values of 1024 random select pixels as x coordinates, 
intensity values of their horizontal neighbors as y coordinates, and intensity values of their vertical 
neighbors as z coordinates. It is obvious that before scrambling adjacent pixels in image Lenna are close 
correlated and thus the plots concentrates on the diagonal direction, namely the line x = y = z. In 
contrast, after scrambling the plot of the same set of adjacent pixels speared almost everywhere. 

n n n 

Results of Ye et al 's |4|,2007, Ye et al 's |7|, 2008 and Abu Dalhoum et al 's |8| are directly taken from listed results in 

^8j. - - 
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Fig. 12. Adjacent pixel correlations before and after Sudoku associated image scrambling for image Lenna; (a) before scrambling; 
and (b) after scrambling. Point colors reflect the distance between a point and the line x = y = z. The closer a point to this 
line, the bluer the color is. 



Adjacent pixel autocorrelation coefficient (APCC) is a common measure used in signal process and 
image encryption [?], |[3|, [ [TT| , [ [T3| , [ [23| , p7| . Mathematically, we can define this autocorrelation 



coefficient as follows: 

p = E[{Xt-^){Xt+i-^)]/a^ (10) 

where Xf denotes the tth element in signal X, E(.) denotes the mathematica expectation as shown in 
Eq. ( [TT] ), fi is the expectation namely the mean value defined in Eq. ([12]), and a denotes the standard 
deviation defined in Eq. ( [T3] ). 

N 
i=l 

fi = E[X] (12) 

a = VE[(X - /x)2] (13) 

Consequently, less correlated adjacent pixels within image X are, the closer p approaches zero. In contrast, 
if adjacent pixels within image X are completely dependent, then |p| = 1. 

Commonly the adjacent pixels in image can be defined along the horizontal direction and the vertical 
direction, respectively. Equivalently, we extract scrambled image pixels along rows and along columns 
respectively and compare the autocorrelation coefficients for these two pixel sequences. Results are given 



in Tables |Vl| and |VII| It is not difficult to see that the proposed method outperforms the three peer 
algorithms again. 
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TABLE VI 

HORIZONTAL ADJACENT PIXEL CORRELATION COEFFICIENTS OF SCRAMBLED IMAGES 



Horizontal APCC 






Test Ima 


ges 






Method 


Cameraman 


Barbara 


Lenna 


#157055 


#69015 


#239096 


Before Scrambling 


0.9565 


0.9238 


0.9709 


0.9522 


0.9613 


0.9796 


Van De Ville et ah 's | 6 |, 2004 


0.6967 


0.6974 


0.7156 


0.6963 


0.7099 


0.7005 


Abu Dalhoum et al 's 2012 


0.0184 


-0.0629 


-0.0204 


-0.0254 


0.1564 


0.1771 


Ye's 125 1, 2010 


0.0827 


0.0398 


0.2130 


0.2270 


0.0359 


-0.0372 


Ours 


-0.0020 


-0.0015 


-0.0021 


0.0024 


-0.0025 


-0.0017 



TABLE VII 

Vertical Adjacent Pixel Correlation Coefficients of Scrambled Images 



Vertical APCC 






Test Ima 


ges 






Method 


Cameraman 


Barbara 


Lenna 


#157055 


#69015 


#239096 


Before Scrambling 


0.9333 


0.8797 


0.9400 


0.9475 


0.9570 


0.9739 


Van De Ville et al 's | 6|, 2004 


0.6613 


0.6699 


0.6741 


0.6690 


0.6785 


0.6865 


Abu Dalhoum et al 's 2012 


0.0543 


0.0172 


0.0476 


0.0376 


0.2162 


0.2314 


Ye's 125 1, 2010 


0.2808 


0.0598 


-0.0955 


0.1015 


0.0860 


-0.0329 


^ Ours 


-0.0033 


-0.0031 


0.0016 


-0.0004 


0.0012 


0.0040 



It is also well-known that whether an observed autocorrelation coefficient p is significantly different 




from zero can be tested by Student's t-distribution p8|-pO|. Specifically in pO[ , it states that "if the 
true correlation between X and Y within the general population is zero, and if the size of sample T, on 
which an observed value p is based is equal to or greater than 6, then the quantity t defined by Eq. ([14]) 

(14) 

is distributed approximately as a Student's t-distribution with T — 2 degrees of freedom where Student's 
t-distribution has the probability density function given by Eq. ( fTS) ), where v is the number of degrees 
of freedom and r(.) is the Gamma function. 



Therefore, given an autocorrelation coefficient of an x i7 image X, it is not difficult to verify whether 
its adjacent pixels are correlated or not by taking the test statistic t for two-sided hypothesis tests. 
Specifically, one can calculate the P-value of a sample correlation coefficient p by finding its P- Value 
under the null hypothesis that the derived statistic i from p follows a Student's t-distribution, where the 
P-value of a given test statistic t is computed as follows 

P-Value(t) = 2 / g{T)dT (16) 
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Statistically speaking, a P- Value is a measure of how much evidence we have under the null hypothesis. 
In other words, the smaller P- Value the more evidence we have against the hypothesis. The range of a 
P- Value is [0,1]. 

In our situation, the sample size T is the number of pixels within image X, which is T = WH. 
Consequently, test statistic t can be derived by the autocorrelation coefficient p and T using Eq. ([14]). 
As a result, corresponding P-value results can be calculated for the autocorrelation coefficients of the 



scrambled image using the proposed method. These results are shown in Table |VIII[ It is clear that 
these P- Values are much greater than 5%, which is a common critical value used in statistics implying 
to reject the null hypothesis if a P- Value is less than 5%, while accept the null hypothesis otherwise. 
Therefore, these results indicate that the correlation coefficient between adjacent pixels in the scrambled 
image using the proposed method are indeed zeros. In other words, the proposed method successfully 
breaks the strong correlations between adjacent pixels after scrambling. 

TABLE VIII 

P-Values of Adjacent Pixel Correlation Coefficients 









Test Images 






Statistics 


Cameraman 


Barbara 


Lenna 


#157055 


#69015 


#239096 


Size 


256 X 256 


256 X 256 


256 X 256 


481 X 321 


321 X 481 


481 X 321 


Degree of Freedom v 


65534 


65534 


65534 


154399 


154399 


154399 


APCC Phorizontal 


-0.0020 


-0.0015 


-0.0021 


0.0024 


-0.0025 


-0.0017 


Statistic thorizontal 


-0.5120 


-0.3840 


-0.5376 


0.9431 


-0.9823 


-0.6680 


Value Phorizontal 


60.84% 


70.10% 


58.08% 


34.56% 


32.60% 


50.42% 


APCC Pvertical 


-0.0033 


-0.0031 


0.0016 


-0.0004 


0.0012 


0.0040 


Statistic tyertical 


-0.8448 


-0.7936 


0.4096 


0.1572 


0.4715 


-1.5718 


F-Ya\uePyerticai 


37.62% 


42.74% 


68.22% 


87.50% 


63.62% 


11.60% 



E. Key Space and Key Sensitivity 

A good image scrambler should have a sufficiently large key space to resist brute-force attacks. Sudoku 
Associated Image Scrambler is designed to be of 192-bit length, which is considered sufficiently large to 
be immune to this type of attacks with respect to the current computer capacities, on the other hand, this 
key space can be easily extended because the number of possible Sudoku matrices are extremely huge. 
A lower bound of the number of 256 x 256 Sudoku matrices is 256! ^ 2^^^^. Considering that majority 
of digital image files are larger than 256 x 256, the number of Sudoku matrices at larger sizes are even 
huge. 

In addition, a good scrambler should have high sensitivity to scrambler key. This has a two-folded 
meanings: a slight key change should lead to significant change in scrambled images during scrambling 
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process, and also significant change in descrambled images during descrambling process. Fig. [13] shows 
the key sensitivity of Sudoku Associated Image Scrambler, with two keys Ka and K^, differentiate from 
each other only for one bit. As can be seen, one bit change in scrambling key leads to two different 
scrambled images, whose difference image is also random-like. And descrambling using an incorrect key 
which is just one bit different from the correct key leads to random-like image. 

Ka = B697F2703EA4347A85D997FB18A1FC3CE7E6901B6A9AE5EA 
Kb = A697F2703EA4347A85D997FB18A1FC3CE7E6901B6A9AE5EA 




VI. Conclusion 

In this article, we mainly discussed the Sudoku associated two dimensional bijections and applications 
for image scrambling. We showed that these bijections are can be defined by Sudoku associated matrix 
element representations, which provide additional and parametric means to denote matrix elements besides 
the conventional way of using the row-column pair. Specifically, these new Sudoku associated matrix 
element representations are row-digit pair, digit-row pair, column-digit pair, digit-column pair, block-digit 
pair, and digit-block pair. 
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Since all these Sudoku associated matrix element representations are parametric with respect to a 
reference Sudoku matrix, it then allows us to denote matrix elements in secret ways and further provides 
Sudoku matrix dependent two dimensional bijections constructed by mapping from one representation to 
the other. We showed that many of these Sudoku associated two dimensional bijections have deterministic 
scrambling effect when we use them for image scrambling. For example, the bijection mapping from 
row-column pair to row-digit pair is equivalent to scramble pixels within a row to different positions that 
none two pixels that originally lies in the same column is still in the same column after scrambling; the 
bijection mapping from block-grid pair to block-digit pair is equivalent to scramble pixels within each 
block and cause a mosaic-like effect. 

Furthermore, we proposed Sudoku Associated Image Scrambler, a simple but effect digital image 
scrambler of using these Sudoku associated two dimensional bijections, by using a scrambling key to 
control these bijections in a parametric way. Because a multiround scrambler is mathematically equivalent 
to a new bijection of composing a series of bijections in each scrambler round, we showed that these 
fundamental Sudoku associated two dimensional bijections can be cascaded together to scramble image 
pixels in a deterministic way. Simulation results of the proposed image scrambler and comparisons to 
recent peer algorithms indicate that Sudoku Associated Image Scrambler outperforms or at least reaches 
state-of-the-art suggested by these peer algorithms ||7|, ||8|, [ |25| with respect to a number of evaluation 



and analysis methods, including human visual inspections, gray degree of scrambling, and autocorrelation 
coefficient of adjacent pixels. Moreover, statistical tests also support that Sudoku Associated Image 
Scrambler does break the strong correlations between adjacent pixels to zero correlations after scrambling. 

Similar scrambling ideas are definitely applicable to other digital formats such as digital audio and 
video. Since data hiding, data encryption, digital watermarking also widely use parametric bijections, we 
believe the proposed Sudoku associated two-dimensional bijections could also be useful in these areas. 
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